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We demonstrate when p-wave pairing occurs between species whose free Fermi surfaces are mis- 
matched the gap generally vanishes over a two-dimensional surface. We present detailed calculations 
of condensation energy, superfluid density (Meissner mass) and specific heat for such states. We also 
consider stability against separation into mixed phases. According to several independent criteria 
that can be checked at weak coupling, the resulting "breached" state appears to be stable over a 
substantial range of parameters. The simple models we consider are homogeneous in position space, 
and break rotation symmetry spontaneously. They should be realizable in cold atom systems. 



I. INTRODUCTION 



Recently there has been considerable interest in a class of possible new states featuring coexistence of superfluid 
O | and normal components. These states arise when there are interactions favoring pairing between fermions that have 
^ . fermi surfaces of different size. For s-wave superfluidity, which has mostly been considered, a breach in the pairing 
occurs in the pairing of fermions which have momenta whose magnitudes lie between two values p±. Two separate 
two-dimensional spherical fermi surfaces, corresponding to gapless modes, open up at \p\ — p±; while paired fermions 
with momenta outside the breach provide a coexisting superfluid condensate. 

The possibility of superfluidity coexisting with gapless states at momenta that span a two-dimensional surface 
originally was suggested by Sarmai^. For spherically symmetric (s-wave) interactions, as he considered, a state of this 
type naturally suggests itself, and a pairing solution can be foundiiSi&iii. The stability of the resulting state against 
phase separation-, or the appearance of a tachyon in the gauge field (negative squared Meissner mass) 7 i 8 i 9 is delicate, 
however. It appears to require some combination of unequal masses, momentum-dependent pairing interactions, and 
long-range neutrality constraintsi2*ii. 

Cold atom systems with mismatched free fermi surfaces naturally arise when the system contains different species 
of atoms, or identical atoms with different spin states in an external magnetic field. Experimental realization of 
fermion superfluids in cold atoms is a major recent development in condensed matter physics^. Manipulation of the 
parameters of Feshbach resonances and various trapping techniques provide control over the form and the strength 
of interaction and effective masses (band structure). In a wider context, having a mixture of bosonic and fermionic 
' atoms and tuning to different parameter regimes, cold atom systems open possibilities to explore new exotic phases^. 
y—( . Most experiments to date have exploited s-wave Feshbach resonances, but p-wave resonance, accessing a state nearly 
t-H ' bound by a dipolar interaction, is also experimentally realizable^. 

At strong enough coupling s-wave breached pairing is realized near Feshbach resonance, opening one or two spherical 
Fermi surfaces^i^, while it is destroyed by instabilities in Meissner mass and number susceptibility^ at weak coupling. 
In this case, the p-wave breached superconductor is a stable ground state as we show in this work. Even when the 
s-wave interaction is repulsive, p-wave pairing instability may arise due to the Kohn-Luttinger effect^, in the same 
way it may support p-wave pairing in liquid 3 He. In QCD with one flavor, the order parameter has total angular 
momentum one, i.e. I = 1 or s = 1, and the preferred phase exhibits color-spin-locking in the relativistic limiljiS. 

Here we suggest analyze a model where the breached pair state arises in the p-wave channel. We find that in 
this context it is quite robustly stable. It seems not unreasonable, intuitively, that expanding an existing (lower- 
ed ■ dimensional) locus of zeros into a two-dimensional surface should be significantly easier than producing, as in s-wave, 
a whole sphere of gapless excitations "from scratch" . In our context, we will show that pairing of mismatched Fermi 
surfaces can expand the locus of gapless states from a line to a torus (polar phase) or from two points into two 
lenticles (planar phase). We shall show that this phenomenon even occurs for arbitrarily small coupling and small 
Fermi surface mismatch, where our mean field approximation should be adequate. We presented a short account of 
some of this work previous!}*^. 

Anisotropic superfluid states that coexist with gapless modes at isolated points or lines in momentum space are 
also well known^ii. Gapless states also are known to occur in the presence of magnetic impurities 21 and, theoretically, 
in states with spontaneous breaking of translation symmetrjiSS, where the gapless states span a two-dimensional 
fermi surface. A similar phenomenon was found in two-dimensional d-wave superconductors subject to an external 
magnetic field 23 . Strong coupling between different bands also may lead to zeros in quasiparticle excitations and 
gapless states^. 

A crucial difference between the model we consider and the conventional p-wave superfluid system, 3 He, lies in the 
dissimilarity of the paired species. Although there are two components, there is no approximate quasispin symmetry, 



and no analogue of the fully gapped B phase 25 arises. 



II. ANISOTROPIC P-WAVE PAIRING 



We consider a model system with the two species of fermions having the same Fermi velocity vf, but different 
Fermi momenta pf ± I /vf- The effective Hamiltonian is 

H = Y}4 a P«p + 4 b -P b -p - A >l b -p + Ap6-pap] (i) 
p 

with = £ p + I, 6^ = £ p — I, £ p = v(p — pf), A p = ^ k Vp_k(a k &l k )- Here the attractive interspecies interaction 
is — Vp_k within the "Debye" energy 2ui£> around the Fermi surface (wd ^ I), and vanishes at larger energies. The 
intraspecies interaction is assumed to be cither repulsive or negligibly small. For the sake of simplicity, we have taken 

the gap function A p to be real. Excitations of the Hamiltonian are gapless, E p = ±^/^p + A 2 + I, provided that 

there are areas on the Fermi surface where I > A p . The gap equation at zero temperature, 

can be simplified by taking the integral over 

A„ . „ / ^V>, „<>A„, ( ta ^ + 6(/ - A n . , ,„ ■ (3) 

Due to hirarchy of scales, wo < Ep, the integration is performed in the narrow window around the Fermi surface 
(BCS approximation), i.e. d 3 k — > vdild^k, where v — l/(2ir) 3 f d 3 k5(£k) = k F /(2ir 2 v) is the density of states. In 
deriving Eq. we neglected dependence of Vp_k on the absolute values of p and k, this is a good approximation as 
long as the "Debye" energy is small compared to the Fermi energy, uj d -c Ep- 

We concentrate on the first harmonic (p-wave pairing) in the expansion of the interaction potential over the spherical 
functions, V(n, n') = g(n ■ n'), with g > 0. The gap equation © allows two solutions corresponding to the value of 
the projection of the angular momentum of the Cooper pair onto some axis z: m = (polar phase), and m = ±1 
(planar phase). We will now analyze these two solutions in detail. 



A. Polar phase: A n ~ Yio(n). 



Let us look for a solution in the form A n = A(z • n) where the direction z reflects a broken rotational symmetry. 
In spherical coordinates, with z || n and 8 = (n, n') the gap equation becomes, 

2 T,„ . „ /Acos0\ , */' „ ; „ fl // + x// 2 -A 2 cosV 



= / d6 sin 6 cos 2 6 In + / dO sin 6 cos 2 6 In . (4) 

vg J V 2uj d J J \ Acos6» J 



where 6* = arccos (I /A), for I < A, and 9* = for I > A. 

The solution of Eq. 10} is different in two cases of small and large Fermi momentum mismatches, 



Small mismatch, I < A. 



Performing the angle integration in Eq. Q), we obtain 
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At zero Fermi momentum mismatch, 7 = 0, the gap is, Ao/2w£> = cxp — -^j w 1.40 exp ^— Using 
dimensionless variables x = J/Aq and y = A/Aq in Eq. (JSJ) we obtain, 



3 

7T / x 



ln- = T - • (6) 



At small the solution of this equation is approximated by y ~ 1 — irx 3 /4. 



2. Large mismatch, I > A. 



The gap equation now takes the form, 



t/2 

■— = / d6»sin6»cos 2 6»l]i 



I + y/P- A 2 cos 2 ( 
2lud 



(7) 



Performing the angle integration we obtain 



± = lm 

vg 3 



I+VI 2 - A 2 
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— — arcsm — 
6 U \I 



(8) 



which in the dimensionless variables x — I/Aq, y — A/Ao is simplified to 
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(9) 



At a; = x c — (l/2)e 1//3 the gap vanishes according to, y w 1.55(x — x,-) 1 / 2 . The solution of the gap equation J^J an d 
@, consists of two branches shown in Fig. 1. 
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FIG. 1: The solutions y(x) of Eqs. 16191 . The lower branch corresponds to the unstable state. Both branches merge at A = I 
(the point B; broken line). Non-zero solutions of the gap equation cease to exist beyond the point A. Stability conditions from 
the energy and the linear response give the points D and D', respectively. Coordinates of the characteristic points (for Ao = 1) 
are A((4/37re) 1/3 = 0.538, e" 1/3 = 0.717), B{ e - n/4 = 0.456, 0.456), (7(0, e~ 1/3 /2 = 0.358), 73(0.475, 0.886), 73' (0.440, 0.917). 



B. Planar phase: A n ~ Yn(n), A n ~ Yi_i(n). 



We now look for a solution corresponding to the momentum m = ±1 of a Cooper pair: A n = A sin (n, z)e ±l< ^, here 
4> is the polar angle in the plane perpendicular to z. The gap equation becomes, 

2 Tj/i • 3/ii /Asinfl\ "f M .3 a , ( I+VP-A 2 sm 2 9\ 

= / d0sin d 01n + / d6 shr 01n — . (10) 

vg J \ 2lo d ) J \ Asm9 ) 

o o v ' 

where 0* — arcsin (I/A), for I < A, and 6* — ir/2, for I > A. Again, we consider separately the regimes of small and 
large Fermi momentum mismatches. 



1. Small mismatch, I < A. 



The first integral in Eq. (ffTTTf) is simple: f^ 2 d9 sin 3 9 In [sin 9] = § In 2 — |. The second integral in Eq. ifTHjl . denoted 
by Ki(I), can be most easily calculated in the following way. Let us first differentiate K\(I) with respect to /, 



arcsin (// A) 

dKi f m sin 3 9 I 1 / I 2 \/A + I 

d9 7W = Wm ~2E? + 4A {' + A~ 2 ) lQ ( A~7 
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Integrating now this equation over / we arrive at, 

i 



(11) 



(12) 



Substituting this expression into Eq. I|10(l we obtain the algebraic gap equation, 

In [1/y] = -*L + £ ( 3 + **) In (I±f) + I In (1 - z 2 ), 



(13) 



where we utilized the same notations as before, y — A/Aq, z = x/y — I /A, x — I / Aq. The gap Ao at zero mismatch 
is Aq/2u>d = 5 exp(5/6 — S/ug) » 1.15 exp (— 3/vg). For x *C 1 the solution to the gap equation has the form, 
y = 1 — 3x 4 /4. Note that the planar phase is more robust than the polar phase with respect to surviving the mismatch, 
in that the gap decreases as the fourth power for a planar phase instead of the third power for the polar phase. 

Qualitatively, the behavior of the planar phase is very similar to Fig. 2 with the following numerical values of the 
characteristic points: xa = 0.674, yA = 0.787, za = xa/va = 0.856; xc = e~ 5 / 6 = 0.435. 



2. Large mismatch, I > A. 



The gap equation at large mismatch becomes 

tt/2 



— = / c20snr01n } = K 2 (I). 

vg J \ 2uj d ) 



(14) 



The integral here is calculated by the method used above [cf. Eq. I|ll|) ]. 



91 J - A 2 sin 2 i 



J 1 / / 2 \, /A + / 
W + 4A 1+ A^ hl T^A 



(15) 



Integrating this expression over /, and noting that K 2 (A) — | In(A/2a;£>) + | In 2 — i|, we obtain, 



i^ 2 (J) = K 2 (A) 



, T dK 2 _ 2 1 / 2A \ 5 J 2 J , 

9/ 3 \2udJ 9 6A 2+ 4Al ' 3A 2 



I 2 \ fA + I 
in 
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(16) 



We observe that the gap equation in the form 



In [1/y] = - 7 -+ 7 -(i + z 2 )\n 



l + z 



1 



-In II - z 2 | 



(17) 



is actually valid for any relation between / and A. We depict the solutions of the gap equation in the polar and planar 
phases in Fig. 2. 



III. STABILITY OF THE OBTAINED SOLUTIONS 



To analyze the stability of the obtained solutions one has to calculate the condensation energy for the two phases. 
The condensation energy in defined as the difference of the thermodynamic potentials in the superfluid Q s and normal 
fL states^i. 




FIG. 2: Solutions y(x) of the gap equation in the polar, m = 0, and planar, m — 1, phases. The lower branch corresponds to 
the unstable state. Non-zero solutions of the gap equation cease to exist beyond the point y'(x) — > oo, where both branches 
merge. The dashed line corresponds to z = 1. 



A. Condensation energy 



The expectation value of the thermodynamic potential at zero temperature is 

tt s = £ p (ul(n Pt+ + rip-) + v%(2 - n p>+ - n p _)) + lA{n p ,+ - n p _) 
p p 

- ^2 V p-k u p v P ( l - n p,+ - n p -)u k Vk(l - n k ,+ - n k _) 
P k 



1 A z 

E o / , p ( : _ n p,+ - "p.-) + E /a K>>+ ~ ( 18 ) 



2 /^ + A 2 



p 



where we distinguish the mismatches I\ , 7 respectively with and without pairing. Here the u p and v p are parameters 
of the Bogoliubov transformation, (it 2 ,, v 2 ) = | ^1 ± Cp/^/^p + i an d n p,± are the Fermi distributions with = 

^ + A p ± Ja ■ The difference in the thermodynamic potentials in the superfluid and normal states is given by, 



n s - n n = l&l - E V^ 2 + A 



which is readily simplified to 



do ri 



n s - n n = v j I - J -f L - - i^il - (A„) 2 e(/ A - |A n |) + i 2 j . (20) 

Note that I — Ia when pairing occurs at fixed chemical potentials, but when pairing is considered at fixed numbers 
of particles the chemical potentials, in general, are different before and after pairing: I Ia- 



B. Polar phase 



Applying the formula (|2U[) to the polar phase at small mismatch, z = I/A < 1, and fixed chemical potentials (so 
I A =7), we obtain 

n s -n n = vA 2 f-l-^ + z 2 ], (2i) 



which is negative provided that z < 0.537. This requires y > 0.886. Thus the upper branch y(x) is stable for 
x = I/Aq < 0.475. This is the analog of Clogston-Chandrasekhar limit, here realized through an anisotropic p-wave 
interaction with the condensate. It is easy to verify that at large mismatches z — I/A > 1, the condensation energy 



n s - n n = vA 1 



arcsm I z 



> 



(22) 



is always positive, meaning that the superfluid state is always unstable. We plot in Fig. 3 the normalized pressure 
difference between the superfluid state with the mismatch / and the normal state with 1 = related to the difference 
at zero mismatch. 



P(I) = 



n fl (/)-n n (Q) 
n s (o)-n„(o)" 



It illustrates that the upper branch is stable for the mismatches less than the critical mismatch, where the dashed line 
crosses the solid line at x = 0.75. The upper branch is unstable for x > 0.75. The lower branch is always unstable. 




FIG. 3: Pressure as a function of the Fermi momenta mismatch x = I/Aq in the polar (solid line) and normal (broken line) 
phases. The upper and lower branches merge at a cusp point (point A at Fig.l). The lower branch corresponds to the unstable 
solution of the gap equation, it is tangent to the normal state at the point where the superfluid gap disappears A = 0. When 
the pressure of the normal state exceeds the pressure of the polar phase (upper branch) , there is a first order phase transition 
at A ^ from the superfluid to the normal state. 



C. Planar phase 



The condensation energy for the planar phase is easily evaluated with the help of the integral, 



arcsm a 



J d9 sin 6 V a 2 - sin 2 



a a 2 — 1 , ( 1 + a 



hi 



I - a 



< a < I, 



(23) 



and yields 



o - <>„ = vA 2 ( - - + — + Z(1 P) In 
1 3 2 4 



1 + z 



1-z 



(24) 



For small mismatches, z < 1 the energy difference is positive for z > 0.623. For large mismatch z > 1, the condensation 
energy is always positive, Q s — f2 n > 0, indicating that the lower branch is always unstable. The dependence of the 
normalized pressure for the planar phase is similar to Fig. 3. 

For / = 0, we obtain the ratio of condensate in the planar phase to the one in polar phase to be e/2 w 1.36, 
indicating that the planar phase is the ground state at zero mismatch. This is in agreement with earlier analyses of 
/ = I pairing 26 ^; the latter contains an extension to higher harmonics. For our specific model Hamiltonian, at weak 
coupling, the planar phase is more stable. 



IV. SUPERFLUID DENSITY, OR "MEISSNER" MASS 



In the s-wave version of our model instability associated with spontaneous breakdown of translation symmetryi& 
can arise. It shows up as a negative superfluid density, or more formally as a negative coefficient multiplying the 



gradient 2 terms in the effective Lagrangian for the superfiuid mode. In the context of superconductivity, the gauge 
invariant form of this term encodes the Meissner mass 2 of the photon. So we can analyze this potential instability by 
checking whether the Meissner mass 2 is positive. 

Following the standard methods in the theory of superconductivity 2 ^ we calculate the super-currents in our system 
under the influence of a homogeneous in space gauge field A, which is assumed to be transverse. The super-current 
is anisotropic, ji — ^^Qik^k with the components obeying Q xx — Q yy . We consider two cases, when the field A is 
perpendicular and parallel to the direction z. It is straightforward to show that the superconducting gap is unaffected 
by the external field in the linear response regime. 

Let us first consider the case A • z = 0. The kernel is given by (for T = 0), 

no on 



2e vipv f f d£du) do n . 9 . , , 

lim ^— / / — sin 2 6 cos 2 6 

q^O to 2 J J 2tt Att 

— oo — oo 

( (iu + £ P+ - I) (ioj + £ P _-I) + A p+ A p _ 



« w + zir + e P+ + ai + )((lo + iiy + e P _ + a 2 _ ) & - / - c P+ ) & 




(25) 



where 8 is the angle between p and z, and <f> is the angle between p and the plane containing the vectors z and A. 
The density of states v = tupf/Sit 2 can be taken independent of / (up to 0(I 2 /E F ) corrections). Since the expression 
(12511 is even in /, the sign of I is irrelevant, so in what follows we take / > for simplicity. The second term in the 
parentheses represents the diamagnetic term (this is easy to verify by calculating first the integral over the frequency) . 
Note, that each term in the parentheses is divergent at large £ and w, but their difference is convergent^. Thus, the 
order of integrals cannot be changed for each term separately, but it can be changed when calculating the difference 
of both terms simultaneously 2 ^. 

It is convenient to take first the integral over <i£ and then the integral over duj, since the second term in the 
parentheses is identically zero when integrating over d£. In the remaining term one can immediately assume q = 
(the integral over <i£ commutes with the transition q — > 0): 

7T OO OO 

— oo — oo 

We will now make use of the following identity (which can be verified by integrating in parts) , 

■ 2 2 2 °° 

^ r^2 , A 2 I /. . , ■ 7^212 = A 2 , 7~. I ■ n2 / d£-p2 , a 9 , 7~. , ■ r\2 ■ (^7) 



\e + A 2 + (w + z/) 2 ] 2 A 2 + (w + il) 2 J + A 2 + (w + il) 



- oo 



The transformation (|27|l improves the convergence at large values of £ and uj. One is now allowed to interchange the 
order of the integrals and calculate first the integral over the frequency with the help of, 

oo 

dbJ 1 



2tt [A 2 + (w + il) 2 } [e + AI + {uj + ilf 




(28) 



Integration over c?£ is now simple and gives, 



OO / v 



A 2 - 



eve+Ai 



= 1- === 6(1-^1). (29) 

v p - a 2 



Substituting this expression into Eq. H26(l we obtain, 



Qxx = ^|^sin 3 flh-_L_ 0(/-|A n |)j (30) 



o 



In the case of A || z, in complete analogy with the derivation of Eq. i|30|) . we obtain, 



(31) 



where the integration over the angle cj) gives 1, instead of 1/2 as for the Q xx . 

Combining both cases, A 1 z and A || z, we can now rewrite the general expression for as follows, 



Szz 

1XX 



e 2 N ( 3 r do n f cos 2 
ml 2 J 47r \ sin 2 9 cos 2 </> J ^/p _ ^2 



e(/-|A n |) 



(32) 



We will now evaluate this expression for the two phases separately. 



A. Polar phase 



The integrals over the angles arc easily calculated to yield, for small mismatches, z = I/A < 1, 



Q, 



Ne 2 



3ttz 3ttz" 



(33) 



The homogeneous superconducting state is stable provided that this expression is positive 1 — 3ttz/A + 3wz 3 /8 > 0. 
This is so as long as z < 0.480, which implies that y > 0.917 or, equivalently, x < 0.440 (see Fig. 4). For larger 
mismatches, x > 0.440, negative values of Q xx probably indicate instability with respect to the transition into the 
Larkin-Ovchinnikov-Fulde-Ferrel state (with paired states of nonzero total momentum). 
For large mismatches, z = I/A > 1, the coefficient 



Qxx — 



Ne 2 



m 



. 3z , / i \ 3z , / 

1 arcsm (z M arcsm z 

2 v ' 4 v 



(34) 



is always negative, indicating instability of the lower branch. 

In the case A || z, for small mismatches z = I/A < 1, the effective density of superconducting fermions in the 
^-direction 



Qz 



Ne 2 



37TZ 3 



(35) 



is positive for z < 0.752, which implies y > 0.717 and x < 0.589. This is a weaker condition compared to the one 
obtained from the density in the x-direction Q xx . For large mismatches A < /, the coefficient 



Ne 2 



1 — - — arcsin (z 1% ) + — \l z 2 — 1 
2 v ; 2 



<0, 



(36) 



again demonstrating that the lower branch is unstable. 



B. Planar phase 



For the planar phase the calculation of the integrals yields for both small and large mismatches, 

l + z 





e 2 N 


y Qxx j 


m 



iT^-|(iT, 2 )m 



1 - z 



(37) 



From these expressions we observe that Q xx reaches zero at z = 0.876, while Q zz always remains positive. Note that 
this critical value exceeds za — 0.856 and zr> = 0.623 meaning that in the energetically favorable state the density of 
superconducting fermions is always positive. 



Qxx 




FIG. 4: Effective density ("Meissner mass") of superconducting fermions in the ^-direction, Q xx , for the polar phase as a 
function of z = I/A. The instability sets in at the point on the upper branch where Q xx becomes negative. The lower branch 
is always unstable. 



V. SPECIFIC HEAT 



A definitive manifestation of the breached states with gapless excitations is the appearance of the term linear in 
temperature in the specific heat, which is characteristic for a normal Fermi liquid. The specific heat is given by the 
formula, 



E, 



dn(E+) t „_dn{E-) 



p dT 



E. 



p dT 



(38) 



where Ep = + ± I, At low temperatures T <C / the first term in Eq. I|38|) gives an exponentially small 
contribution which is negligible. The second term in Eq. 1|38|) provides a leading contribution, 



c- JL f d £ [ d0n ^ 

_oo cosh 



y/e + \A n \*-I 

2T 



(39) 



We will now evaluate this expression for both phases. 



A. Polar phase 



For the polar phase we have, 



cosh 



2T 



(40) 



At low temperatures !"« A the integral over the angle dx can be extended to infinity. Rescaling now xA = x' we 
observe that the integral is conveniently rewritten as the integral over p = \ G, + x' 1 , 



C 



AT 2 A 



dpp- 



(p iy 



cosh 



p-i 

2T 



(41) 



We note that the integrand in this expression is a very sharply peaked function at p ~ I and can, therefore, write, 



dpp- 



(p-iy 



cosh 



p-i 

IT 



I / dp- 



{p-lf IT 3 f dyy 



cosh 



p-i 

IT 



2 J cosh 2 [?//2] 



(42) 



The contribution of the gapless modes to the specific heat is thus, 

n 2 vl 



C 



6A 



-T, 



(43) 



which is a fraction 7/4 A of the specific heat in the normal state. 



B. Planar phase 



The second term in Eq. (|39|l now yields, 



C = 



4T 2 



tv/2 

d£ J dd sin I 





£ 2 + A 2 sin 2 9-1 



cosh 



2T 



(44) 



This integral is a little more tricky for arbitrary ratio between I and A. However, since the principal contribution to 
the integral over dOdS, comes from the area where £ 2 + A 2 sin 2 9 m I 2 , we conclude that for / < A the relevant angles 
are always small 9 <C 1, 



C 



4T 2 



d^ d6l 



£ 2 + A 2 9 2 - I 



cosh 



2T 



(45) 



where again we can extend the limits of the angle integration to infinity due to fast convergence (A ^> T) of the 
integral. Transforming this integral to the polar coordinates similar to the above (p = \/ C 2 + A 2 9 2 ) one obtains, 
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2T 2 A 
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2 (p-if vi 2 f (p-if 2™/ 2 



cosh 
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2T 



2T 2 A 2 



dp- 



cosh 



p-i 

2T 



3 A 2 



T. 



(46) 



This result is valid provided that T « / « A. As might be expected, the linear in T contribution to the specific 
heat is proportional to the area occupied by the gapless modes, i.e. the ~ I/A strip around the equator for the polar 
phase and the islands around the poles for the planar phase. 



VI. COMPETITION OF THE ANISOTROPIC AND MIXED PHASES 



In this section we analyze the stability of the obtained phases with respect to the transition into a spatially 
inhomogeneous state under the condition of fixed number of particles rather than at fixed chemical potentials as 
assumed in the preceding sections. Competition arises between a phase with the finite Fermi momentum mismatch / 
(i.e. with unequal numbers of particles) and a mixed phase consisting of spatially separated regions of superfluid state 
with zero mismatch (equal number of particles) and a normal (unpaired) phase accommodating the extra particles. 
To analyze the stability under the condition of a fixed number of particles one has to evaluate the energy (rather than 
the thermodynamic potential) of both phases. 



A. Mixed phase 



The energy of the mixed phase is given by the sum of energies of the superfluid and normal states, see Ref. (6) for 
the s wave, 

( {^ 2 n)^v (67r 2 n) 2 /3 A 2 \ f (6n 2 n A )^ v {^ 2 n B )^ v \ 
mlx ~ ( '{ 4n 2 2n 2 v 2c ) + X { 8tt 2 + 8vr 2 ) (4?j 

where x/1 — x fraction of the volume is in the normal/superfluid state (thermodynamic limit), total numbers of 
particles (densities) are tia = xTia + (1 — x)n and jib — xt%b + (1 — x)n, the superfluid state is given by the anisotropic 



pairing at zero mismatch containing n particles of each species. We introduced the constant c which assumes values 
c = 3 for the polar phase and c — 3/2 for the planar phase, respectively. There are two variational parameters, x and 
n. It is convenient to introduce the new variables n = ua + Sn and Sn = n B — tia- At small mismatches we obtain 
up to the second order 

{6n 2 n A ) 4 ^v / 2 Sn 1 Sn 2 2 1 - x Sn 2 2 1 - x SnSn\ (6n 2 n A ) 2/3 A 2 , 

E ™ ~ ( v 1+ 3^ + + 9~ nY)~ {l ~ X ' 2tt 2 v ~2^ { ' 

Subtracting the energy of the normal state, 

(67r 2 n A ) i/3 v (6Tr 2 n B ) 4 / 3 v (6Tr 2 n A ) 4/3 v ( 2 Sn I Sn 



and minimizing E m i x — E n with respect to Sn and x, we obtain 



E - ( 6 ^) 2/3 A o f i x * 



2tt 2 u 2c 
(67r 2 n A ) 1/3 «Vc fa 

where A is different for the polar, Ag°', and the planar, Aq', phases. 



,j rn i n ) 



(50) 



B. Energy of the anisotropic phases 

It is convenient to start from a thermodynamic potential at fixed chemical potentials and to perform a Legen- 
dre transformation to find the energy expressed in terms of fixed particle numbers. We begin by calculating the 
thermodynamic potential of the mismatched superffuid state, 

n s (i A )-n = J2 + A p - £ ^ + \ ^ if" A2 ' (51) 

counted from the potential of the normal state at zero mismatch, Qq = £ p — |£ p |; and we use the notation /a 
for the mismatch in the superfluid state. The calculation similar to that from Section III gives, 



Q S (I A ) = n + vj^(-^k- i^ii -A 2 0(1 - |A n |) 



(52) 



We can now perform the Legendre transformation to find the energy of the superfluid state, 

E s = U s H (n A + n B )-\ (n A -n B ) 

= Q +p F (n A +n B )+u J ^ ^-M_/ AV //2 _ A 2^J +lA Sn, (53) 
here Sn = n B — ^a is the difference in the number of particles corresponding to the chemical potential mismatch 7a , 
Sn = ]T e(/ A - V / a + A p ) = 2 I , J ^Ji-A£e(/-|A n |). (54) 

The condensation energy (the difference in the energy of a superfluid state and a normal state) can now be written 
as, 

do n ( A 2 



Es-E n = v J I — y + /a y/1% - A 2 - I j (55) 

where we introduced the chemical potential mismatch / = 8n/2v in the normal state, corresponding to the particle 
number mismatch Sn. Note that in deriving Eq. I|55|) we also used that at non-zero mismatch, E n = Eq + vl 2 and 
ft n — f2o — vl 2 in the normal state. 



C. Polar phase 



Applying the above expression l(55|) for the polar phase, A n = A cos ( 



E K -En = v\- 



6 A 4 



(56) 



where Eq is the energy of the normal state with no mismatch, and the chemical potential mismatch /a is related to 
the particle number difference Sn as, 



2A 



(57) 



We are interested in the case of small mismatches, 4 < Aq, where it is possible to approximate, A/Aq = 1— 7t/ a /4Aq. 
Using these expressions we obtain to the lowest non-vanishing order, 



A 2 

E s - E = 

6 



1-4 



2 / Sn 



7T VfAo 



3/2' 



(58) 



Comparing it with the energy (|50|l of the mixed phase (which is simplified using, v = (67r 2 n y i) 2 / 3 /(27r 2 v) and 



E n 



6 



i - Ve 



Sn 



(59) 



we see that the homogeneous superfluid state is energetically favorable for small mismatches Sn <C ^Aq. 



D. Planar phase 

The energy of the planar phase, A n = Asinfle^^ , counted from the energy of the normal state with zero mismatch 



E s -E = v 



A 2 , l\ AI A 



1 



A 2 



In 



1 + /a/A 



1 - Ia/A 



The particle difference is related to the chemical potential mismatch according to, 



5n = v I 7a - ^ I 1 - ^ ) In 



1 + /a/A 
1 - /a/A 



(60) 



(61) 



For small mismatches Sn <§; ^Ao, one can approximate, Sn = 2^/ A /3A 2 . The gap at small mismatches has the form, 
A/Ao = 1 — 3/ A /4A 4 . Substituting this expression into Eq. I|6l)|l and eliminating Sn with the help of Eq. (|61|l we 
obtain, 



A 2 

E s -E = -v-£ 



1 - 5- 



3 4 / 3 / Sn \ 4/3 



2 7 /3 \vA 0/ 

Again, comparing this expression to the energy H50J1 of the mixed phase (with x m in = V / 3/2(<5n/^Ao)), 



E m ix — Eq = ~ v ~^~ 



1 - V3 



Sn 



(62) 



(63) 



we conclude that the superfluid state is more stable than the mixed state. 



VII. CONCLUSION 



We have presented substantial theoretical evidence that our simple model supports planar phase gapless superfluidity 
in the ground state. For I <C A the gapless modes contribute terms of high order in the mismatch, ~ I 4 in the 
solution of the gap equation and ~ I 2 in the heat capacity, i.e. they represent small perturbations. The planar 
phase is symmetric under simultaneous axial rotation and gauge (i.e., phase) transformation. The residual continuous 
symmetry of this state, and its favorable energy relative to plausible competitors (normal state, polar phase) suggest 
that it is a true ground state in this model. Also, we obtain a positive density of superfluid fermions, suggesting that 
inhomogeneous LOFF phases are disfavored at small /. Direct calculation of the energies in the anisotropic superfluid 
and mixed phases shows that the p-wave breached superfluidity is energetically favorable and the phase separation 
does not occur at small mismatch. 

In some respects the same qualitative behavior we find here in the p wave resembles what arose in s wave^fi. Namely, 
isotropic s-wave superfluidity has two branches of solution: the upper BCS which is stable and - for simple interactions 
- fully gapped, and the lower branch which has gapless modes but is unstable. The striking difference is that in p 
wave the upper branch retains its stability while it develops two-dimensional lenticular surfaces of gapless modes. 
Specifically, the anisotropic p-wave breached pair phase, with coexisting superfluid and normal components, is stable 
already for a wide range of parameters at weak coupling using the simplest (momentum-independent) interaction. 
This bodes well for its future experimental realization. 

In our model, which has no explicit spin degree of freedom, gapless modes occur for either choice of order parameter 
with residual continuous symmetry. By contrast, for 3 He in the B phase the p-wave spin-triplet order parameter is 
a 2 x 2 spin matrix, containing both polar and planar phases components, there are no zeros in the quasiparticle 
energies, and the phenomenology broadly resembles that of a conventional s-wave stated; in the A phase (which 
arises only at T ^ O^S) the separate up and down spin components pair with themselves, in an orbital p wave, and 
no possibility of a mismatch arises. 

Application of anisotropic breached superfluidity to high density QCD may lead to a viable stable phase at low 
chemical potential. Though the gap with higher orbital harmonics is suppressed, the smallest gap defines neutrino 
emission properties and hence cooling rate of the neutron star—. 

It is possible that the emergent fermi gas of gapless excitations develops, as a result of residual interactions, 
secondary condensations. Also, one may consider analogous possibilities for particle-hole, as opposed to particle- 
particle, pairing. In that context, deviations from nesting play the role that fermi surface mismatch plays in the 
particle-particle case. We are actively investigating these issues. 
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